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Abstract 

We present the calculation of the probability production of an 

electron-positron pair in the presence of a strong magnetic field with 

^ ' time-varying strength. The calculation takes into account the pres- 

^ ence of a strong, constant and uniform gravitational field in the same 

direction of the magnetic field. The results show that the presence of 
the gravitational field in general enhances very much the production 
of pairs. In particular, high-energy pairs are more likely produced in 
the presence of the gravitational field than in Minkowski spacetime. 
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1 Introduction 

We have already calculated the production probability of e~-e^ pairs in the 
presence of an overcritical (3> B^t. = rn?c^ /he ~ 4.4 ■ 10^^ G with m and 
— e < the mass and the electric charge of the electron respectively), uniform 
and time varying magnetic field [HEIE] by studying different kinds of time de- 
pendence such as the rotating magnetic field or the fixed-direction magnetic 
field with time-depending strength. Magnetic fields with these characteris- 
tics may be present around highly magnetized neutron stars {magnetars) or 
black holes [H O [6] (obviously, the spatial uniformity of the magnetic field 
is assumed to hold in microscopic length scales of the order of the Compton 
length) where, on the other hand, it is generally believed that gamma-ray 
bursts are produced [the literature on gamma-ray bursts is endless, we quote 
the recent reviews [7]]. Concerning this fact, the production of electrons 
and positrons is in our mind an intermediate step toward the production of 
photons through e~-e~^ annihilations [8] or bremsstrahlung [^, because the 
photons give the only experimental signal that can be observed. It is also 
worth stressing that our approach is always "microscopic" that means that 
we do not consider macroscopic effects like the detailed form of the mag- 
netic field around the astrophysical compact object or the presence of other 
particles or of plasmas in the production region. 

Now, it is clear that in the astrophysical scenario we refer to the gravi- 
tational field produced by the compact object (neutron star or black hole) 
can play an important role. As general aim of our investigation we consider 
the production of particles by a non stationary magnetic field, so we are 
interested in situations where the gravitational effects are not the dominant 
dynamical feature. However, even with these limitations we find interesting 
to consider situations where the gravitational field is strong enough to put 
into the game some aspects of general relativity. Consequently, the simplest 
realistic configuration we may find, i.e. the Schwarzschild metric, was chosen 
as starting point and then, analogously to the magnetic field, the gravita- 
tional field was taken as uniform over a Compton wavelength of the produced 
particle. If one is not too close to the event horizon of the compact object 
the gravitational effects may be treated perturbatively and this was already 
done in pjO] . In order to complete the investigation we suppose here that the 
particle production takes place near the event horizon where, because of the 
singularity of the spacetime metric, a perturbative approach is inapplicable. 
By restricting our attention to the pair production around black holes, the 
investigation is still possible because we can approximate the Schwarzschild 
metric in a form [the Rindler metric [H]], where the general covariant Dirac 
equation is solvable, even in the presence of a uniform magnetic field, pro- 



vided the magnetic and the gravitational fields are parallel 1121 [13] ■ In the 
previous paper [TU] where we have treated perturbatively the gravitational 
field we had the possibility to choose more freely the mutual configuration 
between the two fields. The calculations are performed by using the quantum 
field theory in curved spacetimes [TU [IHl [H] in order to include the presence 
of the background gravitational field and the adiabatic perturbation theory 
[TTj to deal with the time dependent magnetic field. In fact, the presence 
of the gravitational field is taken into account only in the calculation of the 
one-particle electron and positron states without considering the particle pro- 
duction induced by the gravitational field itself as in [181 [13 EDI El] where 
the authors consider the production of particles by a time depending stellar 
[TKl \W\ EH] or cosmological [21] gravitational field. 

The paper is structured as follows. In the next Section, starting from 
the Schwarzschild metric we expand it around the event horizon getting in 
this way a Rindler metric where a particular magnetic field is introduced. 
In Section [3] the Hamiltonian, the one-particle eigenstates of a Dirac parti- 
cle are displayed, with a particular attention to the energy spectrum which 
has, as expected, very different features with respect to the case where only 
the magnetic field is present. Finally, in Section H] the pair production due 
to the time dependence of the magnetic field (in the constant gravitational 
field) is calculated and some conclusions are presented in Section [5l An ap- 
pendix contains mathematical details of some results only quoted in the main 
text. To make easier the reading of the paper some intermediate calculations 
have been omitted: more details can be found in ArXiv under the number 
|liep-ph/0406251. 

Natural units {h = c = 1) are used below and the Minkowski spacetime 
metric is rj^^ = diag(-|-l, —1, —1, —1) with /x, z/ = 0, 1, 2, 3 (Greek indices are 
supposed to run from to 3 while Latin indices from 1 to 3). 

2 Theoretical model 

This Section is devoted to the description of the theoretical model within 
which our calculations are performed. As we have said in the Introduction, 
both the gravitational and the magnetic fields are assumed to be classical 
(not quantized) fields with a given temporal evolution. On the contrary, 
in order to describe the pair creation process, the electron-positron Dirac 
field has to be quantized. Even if we assume that the spacetime structure 
of the gravitational and the magnetic fields is given, it is very difficult to 
determine them because the system built up by the Einstein equations and 
the general covariant Maxwell equations should be solved. To the sake of 



clarity, we remind that we want to calculate the production probability of 
an electron-positron pair in the presence of a static strong gravitational field 
and of a time- varying strong magnetic field. In principle, a possible physical 
system that can give rise to these kinds of fields is represented by a rotating 
black-hole surrounded by a magnetized accretion disHj. In fact, this system is 
believed to be one of the possible candidates as a central engine of the gamma- 
ray bursts and the presence of a strong (~ 10^^ G) magnetic field makes 
possible the energy "extraction" from the black hole through the Blandord- 
Znajek mechanism [23l EU [25l [26]. In particular, the capture of a near star 
by the rotating black hole can lead to the creation of the accretion disc and to 
the production of a strong transient magnetic field |27] [see also the Ref. [28] 
where the time evolution of the magnetic field around a rotating black hole is 
explicitly discussed] . Of course, it is impossible to determine analytically the 
gravitational and the electromagnetic field of this physical system. For this 
reason we will make a number of simplifying assumptions to proceed. First, 
since we imagine the pair to be produced near the black hole event horizon 
we can neglect the gravitational field produced by the accretion disk and 
by the magnetic field and then we can assume that the spacetime metric is 
determined only by the black hole. Even if the black hole is rotating we can 
limit ourselves to study the production of pairs along (or near) the black hole 
rotational axis in order to neglect the effects of the rotation on the spacetime 
metric. In this simplified scenario our starting point is, as in [TU], the metric 
tensor of a non-rotating spherical body of mass M outside the body itself. 
If we indicate with t, X, Y and Z the so-called isotropic coordinates, this 
metric tensor can be written as ESI 



g^y{X,Y,Z) = diag 
where 



(ll.i) " 

F±(X, Y,Z) = 1± "'^ = (II.2) 

^ ^ 4VX2 + Y^ + Z^ 

with re = 2GM the gravitational radius of the body and G the gravitational 
constant. We have chosen the isotropic metric instead of the usual and 
equivalent Schwarzschild metric because from Eq. (III. II) we see that the 



spatial metric is proportional to the Euclidean one and this will simplify 
our future calculations. In particular, if we imagine to pass to spherical 



^Because of the so called "no-hair" theorem [see e. g. [12]] j an isolated black hole cannot 
generate a time-depending magnetic field because its behaviour is completely determined 
only by its mass, charge and angular momentum. In the most general case, a charged 
rotating black hole can generate a static dipole magnetic field. 



coordinates ^29j and we consider two points Pj = {Xj,Yj, Zj) with j = 
that have the same angular coordinates but such that Rj = ^ /Xf + Y^ 



1,2 



and i?2 > Ri then the physical distance A/ between Pi and P2 is given by 
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(11.3) 



with AR = R2 - Ri. 

Now, the e~-e''" pair production is a microscopic process which takes 
place in a volume with typical linear length of the order of the Compton 
length X = 1/m. In this length-scale the gravitational field produced by 
a macroscopic object does not vary very much and the form of the metric 
tensor ( ill. II) can be simplified. As we have said in the Introduction, we want 
to consider here the case in which the pair is created near the black hole 
event horizon lying at \/X^ + Y^ + Z'^ = rc/A q We choose the reference 
system in such a way the pair is created in a volume centered on the z-axis 
then the previous considerations allow us to expand the metric tensor ( III. II) 
around the point Pq = (0,0,rG/4). If P = {x,y,rG/4: + z) with 2; > is a 
generic point near Pq then 
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(11.4) 
(11.5) 



It is clear that we are only interested in the pairs created in the {z > 0)- 
region because those created in the {z < 0)-region will fall into the black 
hole. Concerning the footnote 2 it is also worth stressing that in order that 
our treatment holds the quantity z is not restricted to be microscopic: it is 
important only that z <^ re- 

Now, if we keep only the lowest order non-zero term in g^^{P) then the 
initial metric tensor (III.ip can be written approximately as 



9,u{P) 



9'^{z) 



diag 



2;zY 
re) 



-16,-16,-16 



(11.6) 



^We do not consider the particular case in which the pair is created just on the event 
horizon of the black hole because in this case the direct particle production induced by the 
gravitational field can be the dominating production process and here we are not interested 
in it. Since the process of the direct production of a massive particle can be interpreted as 
a tunnel effect it is sensible to afhrm that it is relevant at distances from the event horizon 
of the order of the Compton wavelength of the created particle. 



This metric tensor has the same form of a Rindler metric tensor describing an 
observer uniformly accelerated in the z direction [TT] |j. Actually, the physical 
meaning of our coordinates is very different from that of the coordinates in 
the Rindler spacetime. For example, while here the coordinate t is precisely 
the time coordinate in the region far from the black hole, the time coordinate 
in the Rindler spacetime is a combination of the Minkowski time coordinate 
and of the Minkowski spatial coordinate along the acceleration. Nevertheless, 
the fact that the two metric tensors have the same form allows us to conclude 
that the metric tensor (111.61) describes a constant and uniform gravitational 
field in the z direction. Observe that no assumption is needed about the 
strength of the gravitational field itself. 

Now, we pass to the mathematical description of the magnetic field. We 
will deal with a uniform magnetic field B(t) with constant direction and 
time-varying strength: 

B(t) = f I . (II.7) 

As we have said in the Introduction, we consider this particular case of a 
magnetic field in the same direction of the gravitational field because only 
in this case the general covariant Dirac equation has been solved explicitly 
[121 [13] and, in order to apply the adiabatic perturbation theory, we need the 
exact solution of that equation [T^. Also, the magnetic field is considered to 
be strong (-B(t) ^ Bcr = m? /e) and slowly varying {B{t)/B{t) <C X~^ = m). 
We choose the electromagnetic gauge in which the vector potential y4^(r, t) 
that gives rise to B(t) is given by 

Ao(r,t) = 0, (II.8) 

^i(r,t) = ^[rxB(t)L = ^, (II.9) 

1 xB(t) 
A,{r,t) = -[rxB{% = ^, (11.10) 

A3(r,t) = i[rxB(t)], = 0. (11.11) 

In the following we will also use the three-dimensional vector A{x,y,t) = 



^We could have scaled the spatial coordinates in order to have exactly a Rindler metric 
tensor, but we prefer to work with x, y and z that are the Cartesian coordinates at infinity. 



{A^{y,t),Ay{x,t),0) with 



Ay[X, t) 



yB{t) 

2 ' 
xB{t) 



(11.12) 
(11.13) 



which is such that d x A{x,y,t) = B(t). 

In order to calculate the pair production probability, we have to build the 
second quantized Hamiltonian of a Dirac field \E'(r,t) in the presence of the 
already introduced gravitational and magnetic fields. We start by writing 
the Lagrangian density of this system which is given by 
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(11.14) 





where g^^\z) = det[glu (z)] is the determinant of the metric tensor, 'y(^^^^(^z) = 
7°ea (z) are the covariant Dirac matrices with Ca (z) the tetrad field 



[HEn] and T\r'{z) = -ia 
tions. The tetrad field Ca ( 



l3u■,^l 



z)e 



'z)/4 are the so-called spin connec- 



AR)o, 



AR)h 
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[z) has been chosen to be diagonal with 
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no sum. 



(11.15) 
(11.16) 



With this choice the spatial spin connections T\ [z) vanish while Fq 



(R), 



Z] IS 



{R) 



independent of z and is given by Fq = 7 7 /4rG. 
The Hamiltonian density is defined as [31] 

^(^) (^, 9,^, ^, 9,^, n, n, r, t) = n(9o^)+(ao^)n-^(^) (^, 9^^, ^, 9^^, r, t) 

(11.17) 
where n(r,t) = 9^(^)/9((9o^) and n(r,t) = (9^(^)/9((9o^) are the fields 
canonically conjugated to \E'(r, t) and ^{v^t) respectively. In our case, it 
can easily be shown that, apart from total derivative terms the Hamiltonian 
density can be written in the form 



^(^) (^, di^, ^, r, t) = 64^"f(r, t)n^^\r, -id, t)^(r, t) 



(H.18) 



where 

Vq I. 4 4 J 4 Tq 

(11.19) 
with P = 7°, ttj. = 7°7^ ay = 7^7^ and a^ = 7^7^ is the one-particle 
Hamiltonian of an electron in the presence of the magnetic field (1II.7P in the 
spacetime with the metric tensor ( ]II.6p . 

If we define, now, the scalar product between two generic spinors ?/'i(r, t) 
and '?/'2(i',i) as 

(^1,^2) = f dS,^-g(^)iz)Mr,th^^^^iz)M^,t) (11.20) 



with S the hyper-surface at constant time then, since ^^/—g^^\z) = 128z/rG, 

(^1,^2) = y"c/r^V^I(r,t)7VgV^2(r,t) = 64y"rfrV^I(r,t)V^2(r,t) 

(11.21) 
and the one-particle Hamiltonian 7^*^^^(r, — iS, t) is Hermitian. This defini- 
tion of the scalar product clarifies the presence of the numerical coefficient 
64 in Eq. (ill. 1811 . We want also to point out here a fact concerning the limits 
of the integrals on the variables x, y and z. In fact, in our model \x\ <^ re, 
\y\ *^ ^G and z <^ tg, but in what follows we will consider only electron 
and positron wave functions that, as |a;|, \y\ or z go to infinity, go to zero 
exponentially with a typical length at most of the order of A then we can 
assume the integrals on x and y as going from —00 to 00 and that on z from 
to cx) without appreciable error. 

Finally, the total Hamiltonian of the system under study is 

Hi'Ht) . / *^""(*. a.*. *. r. .) ^ 64 / .r* V. «)«<«'(. -,0, 0*(r. t) 

(11.22) 
and it depends explicitly on time through the time dependence of the mag- 
netic field. Now, our next step is the determination of the electron and 
positron one-particle instantaneous eigenstates of the one-particle Hamilto- 
nian (llLl9l) . 



3 Determination of the one-particle states 

In this Section we assume that the magnetic field has the same form as in Eq. 
(111.71) but that it does not depend on time. All the quantities, such as the one- 
particle Hamiltonian (III. 191) . that depended on time through the magnetic 

8 



field will be indicated here with the same symbol used in the previous Section 
but, of course, omitting the time-dependence. 

If the magnetic field does not depend on time the eigenvalue equations 

n'^^\r,-id)Uj = WjUj, (III.l) 

n^^^ (r, -id)vj = -WjVj (III.2) 



with Wj, iDj > can be solved exactly [Tzl [13]. Since in these papers a clear 
derivation is presented we will give directly the final form of the eigenstateso 
Now, unlike the Minkowski spacetime case, in the Rindler spacetime it is 
convenient to label the states directly with the energy which is a continuous 
nonnegative quantum number, that will be indicated as E^ independent of the 
others that are the usual nonnegative integers n^ and Ug and the polarization 
cr = ±1. The fact that in the Rindler metric the energy E of the electron 
has continuous eigenvalues from zero to infinity that do not depend on the 
other quantum numbers is the most relevant difference with respect to the 
Minkowski spacetime case [32] • The physical origin of this difference lies 
on the fact that in the present case, due to the presence of the negative 
gravitational potential, the electron mass is no more the lowest energy with 
which the electron or the positron can be created. We will see in Section 4 
how this fact will change also qualitatively the final results with respect to 
the analogous ones in Minkowski spacetime. 

Now, without going into the details we only quote that in the present 
case the electron and positron eigenstates can be written respectively as 



(III.3) 
and 



(T? \ , / fen,rGCOsh(27rErG) 

+ P+Ki/2+2iErGi^kngZ)](pna^ngA^^y) 

(III.4) 

In the previous equations Ku,(x) are the modified Bessel functions with com- 
plex index regular in the limit a; -^ cxo [33] and P± = {1 ± az)/2. Also, the 



''in [T2] another electromagnetic gauge is used but the calculations can be adapted 
straightforwardly to our case. 

9 



quantity 



Kri 



eBn 



m^ 



(III.5) 



can be interpreted as a sort of "transverse" energy of the electron in the 
spacetime with the metric (1II.6P while ipna,ng,a-{x,y) and (^„^ „^ ^.(x,?/) are 
defined as 



/ ^y kn^ + mcxe n,-i,n,{x,y) \ 



^na,n,A^^y) 



2./fe 



v^ 



nd rncr9n^,n,ix,y) 






and as 



'^n^^n.A^^y) 



Zi\ I rCfr) 



^(y \/K, +m aeng,naix, y) 
y/ kng - mae n,-i,nA^, y) 

\ ^kn^+m(T9n„nJyX,y) / 



with the functions 9n^^ng{x,y) given by 



(III.6) 



(III.7) 



leB 1 1 

2TT Udl Ugl 



On„n,{x,y) = jr^_l__l_(at)n,(^t)n,g^p 



eB(x'^ + y'^) 



[111.8) 



and a\ and aj the rising operators relative to the quantum numbers Ud and Ug 
respectively. It can be shown that the states M„^_„g_o-(i?; r) and Vnj^^ng,a{E; r) 
are such that [see Eq. flll.2ip ] 

(Una,ngAE)^^n'^,n'gy{E')) = {Vna,ngAE) , Vn'^^n'g,a'{E')) = 6{E - E')6nay^6ng,n'gSa,a' , 



{Un^,ngAE),Vn'^,n'g,a'{E')) = 0. 



(III.9) 
(III.IO) 



and that the basis {^^^^^^^^(E; r), v„^^„g^^(E; r)} is complete. 

As usual, it is preferable to deal with normalizable wave functions then 
we have to find a convenient boundary condition at a given surface z = h 
that discretizes the energies E. Since the procedure is identical for electron 
and positron states we will consider only the electron states. The functions 
Ki/2±2iErai.'^^na^) go exponentially to zero for large values of k^^z and go to 
infinity as {kn^z)~ for small values of kn^z [33]. For this reason, it is clear 
that 



10 



1. it is impossible to satisfy a "zero" condition for the eigenstates Un^^ng,a{.E; r) 
at a given 4A;„^6 -C 1 or a canonical periodicity condition between two 
points 4:kn^bi <^ 1 and 4,knJ>2 ^ 1; 

2. if we want to build eigenstates with a finite normalization integral 
we have to modify the functions Kii2±2iEra{'^^ndZ) in the region with 

kn^Z < 1. 

In order to do this we proceed as follows. We consider an arbitrary fixed 
surface z = h such that knj^ <^ 1 and assume that the positive-energy eigen- 
states of the one-particle Hamiltonian are the spinors Un^n^.n ,o-(r) defined 
as 



u 



n,na,ng,, 



iv«) 



AtO) 

^ 'n,nd,ng,o 



knjGCOSh{2TTEn,narG) 

47r2 



+P+h/ 



2-2iE„ 



,rG 



{Akn^z)] (Pn^^ngA^^ V) if ^ < & 



kn^rc COs\li2-KEn,nJ'G] 



X 



47r2 

[P.Ki/2+2iE,,,„^rG{^kna 
+P+Kl/2-2iEn,n,rG{4^Ka 






{x, y) ii z > b 
(111.11) 
where Iu]{x) are the modified Bessel functions with complex index regular 
at X = |33], n is a new integer quantum number characterizing the dis- 
crete energies (as we will see these energies will also depend on the quantum 



number Ud) and A^„ 



«) 



na,ng,a- 



and A^, 



(» 



n^nd^Ug^cj 



are two real normalization factors 



to be determined. In fact, by solving step by step the eigenvalue equations 
(IIII.ip one sees that they can be disentangled into an equation depending on 
X and y and another one depending only on z which is indeed a modified 
Bessel equation. The solution space of this last equation is spanned by the 
functions h/2±2iEr^,,^^rG{'^kn^z) and -ft'i/2±2ii<;„,„^rG(4fcnd2;) and we have chosen 
the Kii2±2iEr,,n^rG{'^kn^z) bccause, unlike the /i/2± 2»E„.„^rG(4fc„^^:), they are 
regular in the limit z ^ oo. Now, since Eq. (IIII.ip is a first-order equation 
in the variable z we only require that the spinor M„„^„„n-(r) is continuous 



n^rifi.ng^cj 



at z 
u 



n^rifi^ng.a 



b. By means of this condition and by requiring that the norm of 
r) is the unit, we make the energies discrete and determine the nor- 



malization factors Nn ■, 



«) 



nd,ng,' 



and Nn-i 



(» 



nd,ng,' 



The details of the calculations 
are given in the appendix and here we only quote the final expression of the 
discrete energies [see Eq. (lAlSp 



En,ndrG^Og{knab) 



n— 
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knji 







(III.12) 
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r(<) ov.^ Ar(>) 



and of the coefficients Nn,n^,ng,a and Nn,na,ng,cT [see Eqs. ( 1A5I) and ( lA20p ]: 



8A;„^6 cosh(27rE„,„^rG) Vft^;; 






A^S.,n„. = ^^? = ^ A;n,&-0 (III.14) 



with [see Eq. (IXTOJ) 



2rG 
Qm = log(A;„^6) K^b -^ (III.15) 
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the density of the energy levels. Obviously, all these quantities will be used 
in the calculations but at the end we have to perform the limit 6^0 and 
the physically relevant results must be independent of b. 



4 Calculation of the production probability 

In the framework of the adiabatic perturbation theory, the matrix element of 
the creation process of a pair with the electron in the state uj{r, t) with J = 
{n, rill, ^g, a} and the positron in the state vj'{r, t) with J' = {n', n^, n'^, a'} 
is given by [Hj 

iJi5(t) = {lj{ty,lj'{t)\H(^\tMt)) = 1^^^ [drul{r,t)z{xay-ya.,)vj,{r,t) 

fG J 

where |0(t)) and |lj(t); lj/(t)) = Cj(t)(ij,(t) |0(t)) are the vacuum and the 
pair state at time t respectively [see Eqs. (111.22p . (III.ISP and 0II.19P ]. The 



factor B{t){xay — ya^), using Eqs. (1II.12P and (1II.13P and the expression of 



the induced electric field E(r, t) = —dA{r,t)/dt, can be rewritten in terms 
of the scalar product of the electric field and the velocity operator a that 
is of the work per unit time done by the induced electric field itself. We 
incidentally observe that the electric field E(r, t) is always perpendicular to 
the magnetic field B(t). A more useful form of the previous matrix element 
can be given by using the matrices a± = {a^ ± ioiy)/2 and the relations 



X = \J^^[ag{t) + a\{t) + adit) + 4(t)], (IV.2) 



2V eB{ty 



y = -\i -^Kit) - 4W - «<i(^) + 4W] (IV.3) 



2i y eB{t) 
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with ad{t) [a^it]] and ag{t) [a^(t)] the lowering [rising] operators relative to 
the quantum numbers n^^ and Ug respectively. The result is 

H^'it) = ^^"^^^ [dru^j{r,t)z{a4a,it)+al{t)]-a+[ag{t)+al{t)]}vAr,t). 
rG\/2eB{t) J 

(IV.4) 

Now, in the rest of this Section we will first manipulate the matrix element 

(1IV.4P to put it in the form (lIV.6p . Then we will use it to calculate the 



production probabilities [Eqs. (]IV.14p and (IIV.ISP ] by means of the usual 
adiabatic perturbation theory and finally the total production probability 
per unit volume and unit energy (lIV.34p . 

As we have said at the end of the previous Section, in order to calcu- 
late this matrix elements we should use the expression fllll.lip for uj{r,t) 
with N^fnait) and N^^\t) given by Eqs. (IIII.13D and (IIII.14D respectively 
and the analogous expression for Vj{r,t) U- Actually, an easy power count- 
ing will show that the contribution of the integral on the variable z from 
to 6 goes to in the limit b —>■ 0. In fact, each spinor contains a fac- 
tor l/[kn^{t)b^y\og{kn^{t)b)\ comiug from Ni%{t) [see Eqs. flIII.1311 and 
(IIII.lSp ]. Also, from Eq. ( lA2p we see that the modified Bessel functions 
h/2+2iEn,n (t)rGi4:kn^{t)z) behavc as ^Jkn^{t)z in the integration domain < 
z < b. Finally, because of the presence of the z factor in the matrix element 
( 1IV.4P the result of the integral on z depends on b as kn^{t)b / \og{kn^{t)b) 
and then it goes to zero in the limit 6-^0. In this way, since at the end 
of the calculations the limit 6 — ;> has to be performed, the matrix element 
(1IV.4P can be calculated by using in the whole region z > the expressions 



of the spinors uj{r,t) and vj'{r,t) valid in the region z > b. Actually, we 
can use directly the spinors Un^^ng,a{E;r) and f„/ „/ ,j/(i?'; r) multiplied by 

Nn^'{t) and A^„/ (t) respectively because the presence of the factor z in the 
matrix element flIV.4p makes finite the resulting integral from to oo [see 
also the general formula flA12p ]: 



nrf,no,CT;n' n' o-'V ' ' '' 



na,ng,a;n'^ 



a' 



rGJ2eB{t)gn^{t)gn'^{t) 
X / c^ru],^,„^_,(E;r,t)z{a_[ad(t) + aJ(t)]- 

(IV.5) 



^Remind that these quantities are now time-depending because the magnetic field de- 
pends on time. 
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At this point we have to substitute Eqs. fllll.Sp and flIII.4p with the time 



dependent magnetic field in place of u'l^^^ ,cr(-^5 ^^ ^) ^^^ ^"' ,«' ,o-'(-^'5 ^5 ^) ^^^1 
apply the various operators. The calculations are involved iDut straightfor- 
ward. We will not report them here but we give the final result 



z^(«) (E E'-t) = ^^^(^^ / cosh(27rErc)cosh(27rEVG) ^^ 



^ ^ ^ ^ (^g + 1) [kuM - "^^] [kn,+i{t) - ma'] ^ 

X Re ((1 - ia){l + ia')Ina,na+i{E, E'; t)) (5„^,„;-i5„^+i,„;j+ 



I rid [kn^t) - ma] [kn^{t) - ma'] ^^ 

X Re ((1 - ia){l + ia')^n„n,(^, E'; t)) 5n„n'^5n^,n'- 



'ug [kn^t] + ma] [kn^-iit) +ma^^ 



X Re ((1 + ia){l - ia')In^^na-i{E, E'; t)) 5„^-i,„' (5„^_i,„' - 



'nd [kna{t) + ma] [kngjt ) + ma'] 

X 



X Re ((1 + ia){l - ia')In^,nAE, E'- 1)) K^^n'^K^n^ 

(IV.6) 
where the adimensional function 

^hit) \ ^^ ( Akiit) 



Xl^l'(E,E';t) = / dssKi/2-2iErG , , ^ Ki/2+2iE 

Jo \ \/2eB(t) 



re I ' —"5 



^2^^ J KV^^W) 

(IV.7) 
has been introduced. 

Before continuing we want to point out that from Eq. ( JIV.6P it can be seen 



that the total angular momentum of the electron-positron field is conserved 
in the transition in fact, in any case 

nd-rig- - + n'd-ng + - = nd-ng + nd-ng = 0. (IV.8) 

Of course, this selection rule is a consequence of the fact that the time evo- 
lution of the magnetic field does not break the rotational symmetry of the 
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system around the z axis or, in other words, of the fact that the total angular 
momentum along z and Ti.^^^r, —id,t) commute. 

We also observe that, as in the Minkowski spacetime, it is impossible in 
the case under study to create a pair in which the electron is in a (n^ = 0, cr = 
— l)-state and the positron in a {ug = 0,a = +l)-state [U |3]. In fact, we 
remind that this selection rule holds in general for the eigenstates of o"^ when 
(Tz anticommute with the time-derivative of the one-particle Hamiltonian [3] 
and it can easily be shown that this is true in our present case because the 
gravitational field changes only the longitudinal structure of the one-particle 
electron and positron wave functions. 

Now, since we are interested only in the strong magnetic field regime in 
which eB{t) ^ m^, we can simplify the expression of the transition matrix 
element flIV.6p by taking into account only those transitions whose proba- 
bilities are proportional to the lowest power of iri^ /eB{t). In the framework 
of the adiabatic perturbation theory the first-order transition amplitude in 
B{t) of the creation of a pair at time t in the state with quantum numbers 
{E, Hd, Hg, cr; E\ n'^, n^, a'} is given by [H] 

7n,,n,,.;n'„n',AE, E' ; t) = ^^ ^ dt' H^f^^^^^^^.^^^^^^.^E ^ E' ' t') cxp {^{E + Er) 

(IV.9) 
and the corresponding probability is the square modulus of this number. It is 
evident that, since the energies E do not depend on -B(t), we can perform the 
(m^/ei? (t))-power counting directly on the matrix element (1IV.6I) . To this 
end we need the general behaviour of two particular classes of the integral 
fHV:?!) that is Io,n'^iE, E'; t) with n'^ > and Ina,nAE, E'; t) with Ud > 0. By 
reminding the expression fllll.Sp with the time-dependent magnetic field for 
kn^{t) and by using the general formula flA12p it can easily be seen that 

1/2 

if n;^ > and eB{t) > m^, (IV. 10) 



^o,n'iE,E';t) 



V "'d 



m 



V^m 



In„n,{E, E'; t) iind>0 and eB{t) > m^ (IV. 11) 

nd 

where, for later convenience, we have also pointed out the dependence on 

the quantum numbers Ud and n'^. Obviously, the integrals Xn^,nd±i{E,E';t) 

behave as the integral ( lIV.lip and then this criterion allows us to neglect 



the transitions in which the electron is in a (n^ = 0, a = — l)-state or the 
positron in a {ug = 0, cr = +l)-state [see Eq. (]IV.6p 



Another criterion we will use to select only the most probable transitions 
is the dependence of the corresponding probabilities on the quantum num- 
bers Hd and Ug. As previously, we can work directly on the matrix element 
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flIV.6p by keeping in mind that at the end we will sum the probabilities with 
different values of ria and n^. Now, we have explained in [3] that the internal 
consistency of the model requires that the sum on Ug (and on n^) cannot be 
extended up to infinity but that they must be stopped up to a certain Nuif) 
corresponding through the relation [see also [3j] 



iVM(«) . '-^R^ 



32 



±M 



(IV. 12) 



to a fixed Rj^m whose physical meaning is explained below (the presence of 
the factor 1/32 is due to the fact that the spatial metric in Eq. ( 111.61) is 
only proportional to the Euclidean one). Now, since we are considering the 
production of electrons (positrons) in states also with Ud {rig) different from 
zero, in general the transverse motion of a classical electron (positron) is con- 
fined within a circle with radius 2R±m, then this quantity can be assumed as 
the radius of the quantization cylinder. Coming back to the matrix element 
(1IV.6I1 . we see that it contains two kinds of terms, the first one being propor- 
tional essentially to y/n^ and the second one to ^/nd. By taking into account 
Eq. OIV.lip it is easy to see that the first kind of terms gives rise to final 
probabilities proportional to Nf^{t) log NM{t), while the second one to final 
probabilities proportional to Nlj{t). For all these reasons we can consider 
only the transitions to states with n^ > and n' > and approximate the 
matrix element (1IV.6I) as 



r(-f?.) 



H III (E, E ; t] 

nd,ng,(T;n'g,n^,u'\ ' ' / 



ii3{t) 



(27r)25(t) 



^cosh(27rErG) cosh(27rE'rG) x 



X 



{ug + l)^Jnd{nd + 1) 



QnA^)Qnd+l{t) 
X Re ((1 - i(j){l + i(T')In^^ni+l{E, E')) (5„^,„' -I'^n^+l.n' + 



rig^JudJUd- 1) 

QnAt)Qnd-l{t) 

X Re ((1 + ia){l - ia')Xn^^nd-i{E, E')) 5„^-i,n^5„^_i,„;^ 

(IV.13) 

where we have pointed out that in the strong magnetic field regime if n^^ > 
and n^ > the integrals Xn^^n^±i{E , E') do not depend on time [see Eq. 

(HH]. 

At this point, by inserting this matrix element in Eq. (1IV.9I) . by squar- 
ing and summing on the polarization variable a we obtain the differential 
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probabilities 

dPna,ng;na+i,ng+i{,E,E';t) = Y^i^ig + 1) A/n^^(n7+^ cosh(27rErG) cosh(27rEVG) : 

X '^7/^^f;)f^'' IHE, E'- 1) f dEdE', 



dR 



n^+l,ng+l 



{IV. U) 
.^^^n,{E,E';t) = —{ug + l)^/ndind + I) cosh{27r Era) cosh{27cE'rG) 

\In,+l,n, {E, E') I 1^^^^ ^,_ ^^ |2 ^^^^, 



{E + E'Y 



(IV. 15) 



where we have muhiphed by the number of electronic states Qn^{t)dE with 
energies between E and E + dE and by the number of positronic states 
Qn^+i{t)dE' with energies between E' and E' + dE' and where 

T{E,E'-t) ^ |yt'^^exp{^{E + E')t'). (IV.16) 



We point out that, as expected, the probabilities (1IV.14P and (IIV.ISP do not 
depend on the unphysical parameter b. 

Now, we want to calculate the probability dP{E, E'; t) that a pair is 
present at time t with the electron with energy between E and E + dE and 
the positron with energy between E' and E' + dE'. To do this we have to sum 
on the remaining quantum numbers n^ and Ug. As we already know, both 
the series on n^ and Ug are diverging then we can perform the summations 
by assuming n^ ~ n^ + 1 and rid — Ud + 1 because the most relevant terms 
are those with n^, :» 1 and Ud ':^ I- Starting from Eqs. ( JIV.14P and ( lIV.lSp 
we have 



dPiE,E';t)c^- 



TT" 



NMit) 

E 

"0=1 



n„ 



cosh(27r_ErG) cosh(27r_E'VG) x 



X 



j:Zfnd\Ir.,,nAE,E')f 



{E + E') 



\T{E,E';t)fdEdE' 
(IV. 17) 
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where A^m(^) has been defined in Eq. 0IV.12p . 

The next step is the exphcit calculation of the functions T„^ „^(ii^, E') and 
J-'{E, E'; t) defined in Eqs. ( 1IV.7I) and ( 1IV.16I) . By using the general formula 
(1A12I) and the properties of the F function [33j it can easily be shown that 

l_^^E-E')rG 7r{E + E')rG 



^na,na{E,E') 



n 



(IV. 18) 



2nd cosh {n{E - E')^g) sinh {n{E + E'Yg) ' 

In order to evaluate the function J^{E, E'; t) we have to assign the time 
dependence of the magnetic field. We assume that 



B{t) = Bf + {B, 



B/)exp( -- 



(IV. 19) 



with Bi < Bf and it results 



T{E, E'; t) = / ds' 
Jo 



Bf -Bi 



— e^^{-s' + i{E + E')Ts') 



Bf + {Bi - Bf) exp [-s\ 

(IV.20) 
with s' = t'/r. Now, from a physical point of view we are interested only 
in the cases such that Et ^ 1 and E't ^ 1. In fact, the eigenvalues E 
correspond to the classical energies [29] 

Ei-\z) = -p=^- (IV.21) 

with v^ the square of physical velocity as measured by the local observer in 
the gravitational field. Also, r is a macroscopic time parameter connected to 
the typical evolution time of the black hole. Analytic estimates of the time 
duration of the formation of a black hole suggest that t ^tq |29]. In this 
hypotheses, even if z ~ A then E^^^\z)t ^ 1 because we are interested in 
energetic electrons with a Lorentz factor 1/Vl — v"^ ^ 1- For this reason we 
can give the following asymptotic estimate of the integral flIV.20p 

t/T + i{E + E')t) 1" 

~~B'i 



HE,E';t) 
or, as t 



Bf — Bi 



i{E + E')T 



exp 



Bf + {Bi - Bf) exp {-t/r) 



(IV.22) 



CXD, 



J^{E,E';t^ oo) 



Bf - B, 



(IV.23) 



Bi {E + E')t 
By substituting this expression and Eq. (IIV.ISP in Eq. (jIV.lTp we can write 



the asymptotic value of the probability dP{E, E'; t —>■ oo) as 



dP{E,E';t ^oo) 



eBfR^ 



64 



log 



e5fi?2 



32 



Bf — Bj 

BiT 



r%G{rG]E,E')dEdE' 
(IV. 24) 
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where we have made the substitutions [see Eq. flIV.12p ] 



NMit^oo) 1 / n n2 \ 2 



E n, -^iNUt -.00) = I rJ^] , (IV.25) 



and where the function 

1 + [{E- E')rG? cosh{27TErG) cosh{27rE'rG) 



G{rG;E,E') 



[{E + E')rG]'^ cosh2(7r(^ - E')^g) sinh2(7r(E + E')rG) 

(IV.27) 
has been introduced. We have pointed out the dependence of ^(r^; E, E') on 
the parameter tg because, as we have said, we are interested in energies E 
and E' such that Etg ^ 1 and E'tg ^ 1. In this energy region the function 
Q{rG', E,E') strongly depends even on small changes of E and E' through 
the hyperbolic functions. In particular, it can be seen that 



[{E + E'jrcl^ 



1 



sinh^(7r(E + E')^g) cosh^(7r(^ - ^')^g) 

(IV.28) 



From this expression and by reminding that E,E' > 0, one can easily show 
that 

lim ;;i'f/^'!;°'l g(rc; e, e) ^h(E- e) (iv.29) 

rc^oo 1 + y[E — E'jrGi tt 

and then that 

g{rG; E, E') ~ l J^^ ' ^') if e,E'>Q and Er^, E'tg » 1. 
2-11 [EtgY Tg 

(IV.30) 

With this result and by integrating Eq. (II V. 241) with respect to the positron 
energy E' we, finally, obtain the probability that an electron is created at 
time t -^ 00 with an energy between E and E + dE such that Et > Etg ^ 1 
in the form 

1 (<'B,RlM\\^JeB,Ri„\ /B,-Bi\\adE 



,P(E; < ^ 00) ^ - ^-^^) log ^^^J ^^^ j ^. 

(IV.31) 
In order to obtain a probability per unit volume we have to give an esti- 
mate of the height of the quantization cylinder. Now, we have said that the 
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modified Bessel functions Ki/2±2iErGi^^nd,it)z) (we refer to tlie electron wave 
functions but an identical conclusion can be drawn for the positron ones) are 
exponentially decreasing as kn^{t)z ^ 1. Now, by using the well-known oscil- 
lation theorems it can be shown that if Era ^ 1 the exponential behaviour 
of the function Ki/2±2iErGi^^nd(t ~^ 00)2) starts at 

2ErG 2ErG 

zq ~ ^, , ^ — - < , ^ . 1V.32 

For this reason we can assume 

as the volume of the quantization cylinder and then 

dP{E;t^ 00) I feBA "'^ ( eBfRJ.A f Bj - BA ' ^ 

dVdE TT V 128 y ^ V 32 J \ rBi J E^ ' 

(1V.34) 
We first observe that the disappearance of the electron mass and of the 
gravitational radius of the black hole is due only to the fact that we are 
working in the strong magnetic field regime and in the high-energy region. 
In order to clarify the meaning of some variables appearing in Eq. flIV.34p it 



may be useful to recall the physical model which is the starting point of the 
whole investigation. There is a bundle of lines of magnetic flux that is limited 
transversally but its size is very large at the microscopic scale where the field 
may be considered uniform. We remark that the existence of two scales, one 
of astrophysical origin and one set by the elementary particles is essential 
in the whole treatment. The time variation, slow at microscopic scale, may 
involve both changes in strength and in direction of the magnetic field and in 
the present case, as in |1J, only the change in strength is considered. In any 
case the variation of B(t) gives rise to an electric field which is not uniform 
in the transverse variables: starting from the center of the bundle the electric 
field produced by the magnetic- flux variation is [r x B(t)]/2, therefore the 
factor {Bf — Bi)R^M/'T gives the order of magnitude of the induced electric 
field at the boundary of the volume. Concerning this fact, we point out that a 
meaningful interpretation of R^m can be given. In fact, R_lm can be seen as 
the typical length scale in which the magnetic (gravitational) field produced 
by the astrophysical compact object can be assumed to be uniform. 

Also, before comparing Eq. (1IV.34P to the analogous result in [T] where no 
gravitational field effects were taken into account we have to stress that the 
comparison can be only qualitative because in [Ij a different kind of magnetic 
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field time-dependence was used. Nevertheless, we note in the present case a 
qualitatively completely different dependence on the physical quantities like 
the magnetic field strength or the electron (positron) energy. In fact, here 
the probability per unit volume depends on the (5/2)-power of the magnetic 
field strength while in [1] it depended on the (3/2)-power of the magnetic 
field strength. Also, the probability production per unit energy scales here 
as E~^ in the high-energy region while in the preceding case the probability 
behaved as E~^. In this way, a general enhancement of the pair production 
is observed and the production of high-energy pairs is strongly favored in the 
presence of the gravitational field then the effects of the gravitational field in 
the pair production process are really relevant and they can not be neglected. 
To give a quantitative estimate we first integrate Eq. flIV.34p from E„ 



100 r^^ (remind that we assumed Ere ^ 1) to infinity. By indicating the 
resulting total probability per unit volume as dP(t -^ oo)/dV, we obtain 



dV TT V 64 y ° V 32 J \ TBi J El 

Now, by assuming re = 4.4 Km as for a 1.5 solar masses black hole and r = 
10~^ s > re, it also results EmT ^ 1. By using the typical values R^_m = ^Gi 
Bi = Bf/10 = 10^^ G and by assuming the magnetic field parameters used in 
m to be such that Bq = Bi, B{t) ~ Bf and b = {Bf — Bi)/T we obtain that 
the ratio between Eq. (20) in P (divided by the volume ZnR^j^) and the 
previous probability is of the order of 10~^^ that clearly shows how relevant 
are the effects of the presence of the gravitational field in the pair production 
process. 

Finally, by using the formula 



dVdE TT V 128 y ° V 32 J \ rBi J E, 

(IV.36) 
we give an estimate of the total energy produced through the mechanism 
at hand. By substituting the previous numerical values and V^^^ = r^ we 
obtain S ~ 10^^ erg. This result must be considered with some care because 
it is roughly two orders of magnitude larger than the rest energy of a black 
hole with a mass 1.5 times the solar mass (on the other hand, the magnetic 
energy associated to the magnetic field is negligible with respect to the black 
hole rest energy). One reason of this overestimate is that for simplicity 
we have assumed the pair production to be isotropic in space and we have 
obtained it by integrating on all the possible directions. Instead, actually, 
the strong magnetic field deflects the electrons and the positrons in such a 
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way their "macroscopic" motion is along the magnetic field itself (of course, 
the particles also rotate in the plane perpendicular to the magnetic field but 
along circles with a microscopic radius) . Another reason of the overestimate 
is the fact that we have not considered here the backreaction of the created 
pairs on the existing magnetic field. This backreaction reduces the magnetic 
field strength (and then the final pair production yield) because the created 
electrons and positrons rotate in such a way as to produce a magnetic field in 
the opposite direction of the strong magnetic field B{t). It must also be said 
that the value £ ~ 10^^ erg would be too large with respect to the typical 
energy carried by a gamma-ray burst which is ~ 10^^ — 10^^ erg. In fact, as 
it has been shown in |35], at so high luminosities as those considered here 
~ S/t ~ 10^^ erg/s, almost all the pairs produced annihilate into photons. 
Concerning this fact, we have studied the spectrum of the photons produced 
in Minkowski spacetime as a consequence of the annihilation of the pairs 
created [H] or as synchrotron radiation [9]. Analogously to what we have 
seen in [8], we expect that the annihilation spectrum is again peaked around 
the electron mass but less sharply because the creation of high-energy pairs 
is here less suppressed than as in Minkowski spacetime. Also, we have seen in 
[9] that the spectrum of the photons emitted as synchrotron radiation by the 
pairs already created has some qualitative similarities to those of gamma-ray 
bursts but that it decreases at higher energies faster than u~^ which is the 
typical behaviour of the high-energy part of the gamma-ray bursts spectra 
[7] . In the present case we expect qualitatively a less fast decreasing because 
the high-energy pair production probability scales as E~^ instead that as 
E~^ as in Minkowski spacetime. 

5 Conclusions and final considerations 

In this paper we have calculated the probability that an electron-positron 
pair is created in the presence of a strong magnetic field with time-varying 
strength and of a constant and uniform gravitational field parallel to the 
magnetic field and represented by the Rindler spacetime metric. We have 
compared this probability with the analogous one calculated in Minkowski 
spacetime and we have shown how large are the effects of the gravitational 
field. From a physical point of view this can be expected because the pair is 
created here in the presence of a negative gravitational potential in such a way 
the rest mass of the electron does no more represent an energetic lower limit. 
Concerning this fact, we observe that the presence of the gravitational field 
makes possible the creation of e~-e~^ pairs that cannot fiy to infinity because 
they do not have enough energy. What we really expect is that the charged 
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particles created with such energies annihilate inside the gravitational field 
giving rise to photons which may fiy away. So, in the presence of the gravita- 
tional field, the spectrum of photons produced through annihilation extends 
to low frequencies. Actually, this spectrum is also expected qualitatively to 
be broadened at higher frequencies with respect to the Minkowski spacetime 
case because high-energy pairs are more easily created in the presence of 
the gravitational field. In fact, being E the electron (positron) energy here 
the production probability per unit energy scales as E^^ while in Minkowski 
spacetime as E'''^. 

We have also pointed out a qualitative different dependence of the pro- 
duction probability on the magnetic field strength B: here it is proportional 
to S^/^ while in Minkowski spacetime to B^^^. In this respect we have ob- 
served a huge enhancement in the production probability with respect to the 
Minkowski spacetime case. Concerning this effect of the gravitational field, 
it is worth considering the "intermediate" situation in which a gravitational 
field is present in the production region but it is not so strong as in the vicin- 
ity of the event horizon of a black hole. This situation happens, for example, 
if the pair is imagined to be produced around a neutron star. We have al- 
ready considered this case in [10] and we refer the reader to that reference 
for details. We have seen there that, in the case of a neutron stars, the ef- 
fects of the gravitational field can be treated perturbatively. In the different 
configuration in which the magnetic field was rotating we have shown that 
the presence of the gravitational field in that case enhances the production 
of pairs but does not change the scaling of the production probability with 
respect to the magnetic field. This conclusion is not unexpected in a pertur- 
bative regime and it is sensible to affirm that it also holds for the present case 
of magnetic field varying only in strength. In this respect the conclusion is 
that the strong enhancement in the pair production yield we have obtained 
here is due to the non-perturbative effects of the gravitational field and it 
only occurs for the case of objects with an event horizon. Correspondingly, 
by means of a numerical estimate we have shown that the order of magnitude 
of the energy produced through this mechanism is very large and (at most) 
~ 10^^ erg in the physical conditions previously discussed. 

We want to conclude by comparing qualitatively the pair production 
mechanism discussed in the paper with the direct pair production by a grav- 
itational field. In fact, as we have said in the Introduction, many works have 
been done about the production of particles (photons, e~-e+ pairs) in the 
presence of a time-depending gravitational field [201 [2IJj- We quote, in par- 



^Actually, the creation of particles by a strong gravitational field is not bounded to the 
case of collapsing objects, but it can also happens for static situations (see, for example, 
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ticular, the seminal paper [19] by Hawking where it was shown for the first 
time that a collapsing black hole with mass M emits electrons (positrons) 
with a thermal energy spectrum proportional to the Fermi-Dirac factor 

^ (V.l) 



ex^iuj /ksT) + 1 



where uj is the energy of the electron and ksT = 1/{8ttMG) = l/(47rrG') 
with ks the Boltzmann constant. From the previous equation it is clear 
that the main emission of particles takes place for wavelengths of the order 
of the Schwarzschild radius of the emitting object. Now, as it has been 
already stated, we have studied the creation of particles with microscopic 
wavelengths of the order of A = 1/m, so the process under study takes places 
in a region of energies where the emission for pure gravitational effects is very 
small for every astrophysical object. The situation could become different 
when the particle production is primed by a rotating black hole. In this 
case, the quantization of a field in the corresponding Kerr metric is a very 
complicated issue because the definition itself of a vacuum is problematic 
[36] . Nevertheless, also in [19], it is argued that the energy spectrum of the 
particle created is the same as in Schwarzschild case [see Eq. fIV.ip ] but with 
the particle energy u substituted by cj — miVtK where mi is the quantum 
number of the particle angular momentum component along the black hole 
rotational axis and VLk is the angular frequency of the black hole. In this way, 
we can conclude quahtatively that the possibility of coexistence of the two 
emission mechanisms could be found only in the case of really huge values of 
the angular momentum of the emitted particle but a quantitative estimate 
would involve detailed and non trivial calculations. 

Appendix 



In this appendix we will impose that the spinor (jIII.lip is continuous aX z = h 



and that its norm is unit. As a result, we will discretize the energies E and 
determine the two factors Nn,na,ng,a and Nn,nj^,ng,a appearing in Eq. (1III.11I) . 
The continuity condition is satisfied if 

Ni':lh/2+2^E„,„^rai'^kn,b) = N>lK„^+2^E^^^^ra{^knJ)) (Al) 

where we pointed out that NIC}1^ and A''^^^^ can not depend on Ug and a 
and that the energies depend on a new integer quantum number n. Since 



the production of particles around the so-caUed "eternal black holes" in [14]). 
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kna^ <^ 1 we can use the approximated expressions of the modified Bessel 
functions near the origin [33j 

/l/2+2i£;„,„.rG(4fc„/) 



1 



r(3/2 + 2tEn,n,rG) 



i^Kb) 



l/2+2iE„,„ re 



and Eq. (JAll) becomes 



(A2) 
(A3) 



l—AiEro 



(A4) 



where the property T{z+1) = zT{z) has been used. By equating the modulus 
and the phase of the left and right hand sides of Eq. flA4p we obtain the two 
real conditions 



'^'"^ 8K,b cosh(27rE„,„,rG) "•'^''' 

1 



(A5) 



arctan (4E„,„^rG) + 2 arg ( T ( - + 2iE„,„^rG ) ) - 4E„,„^rG log(2A;„^6) = 2mT n = 0, ±1, 



(A6) 



where the following property of the F function has been used J33] : 

2 



r|i + >« r i->« 



r|5 + >? 



vr 



cosh(7r^) 



with ^eR. (A7) 



The condition (1A6P determines the allowed discrete energies while, in order 
to determine NnJ^^ "^^ have to require the following normalization condition 
[see the expression (III.2ip of the scalar product] : 



64 / rfr4,„^„^ ,^(r)u„,„,,„^,,(r) = 1. 



(A^ 



It can be seen that the previous condition is equivalent to require that 
/cn^rc cosh(2£;„,„^rG) '" '"^ 



64- 



47r2 



Jo 

dz\Ky2+2^Er.,^^rai^K,z)\' 



(A9) 



1. 
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By using the approximated expressions (IA2p calculated in Akn^z the first 
integral gives 



Jo 



|2 1 cosh(27rE„,^^rG) {K^by 

A;„^ 1 + (4E„,„^rG)2 re 



(AlO) 



The second integral can be evaluated by using the following identity 

1 r /'oo 

dz |i^l/2+2i£;„,„,rG(4fc„,2;)|^ = lini / dz{4:kn^zy \Ky2+2iE„,„^rGi4:KaZ)\'^ - 

^~*^ VJo 

- / dz{Akn^Zy\Ki/2+2iE^,„^rG{^kn^z)\ 

Jo 

(All) 

The first integral in the right hand side of this equation is a particular case 
of the general formula [31] 







r(i-p) 
X r ( i^^±A^^ r (l^I^A^^ X 



l-p + A + /il-p-A + /i^ ^ 6^ 

2 ' 2 ''-^''-^ 

(A12) 



where Re(a + 6) > and Re(p) < 1 — |Re(A)| — |Re(yu)| and where F(r, s; m; ^) 
is the hypergeometric function. Instead, in the second integral in the right 
hand side of Eq. ( JAllj) the approximated expression ( 1A3|) of the modified 
Bessel function calculated in 4fc„^2; can be used. Obviously, even if these 
integrals are both diverging in the limit e — > 0, their divergences must cancel 
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each other because the left hand side of Eq. flAlip is finite. In fact 



hm 



" ^oo j'b 

j dz{Akn^zY\Ki/2+2iEr..r.^rGi.'^K^z)\^ - / dz{AK^zf \Kil2+2iEr.^^^rGi.'^K^z)\^ 

Jo Jo 



hm 



2-2+^ ^/2 + e 



r(i + 6) 



r ( ^ + 2tE„n,rG 



r 



1 

Ye 



r I - + 2iEn.ajG 



(4A:„,6) 



Ak 



nd 



r ( - + 2iE„„ndrG 



hm 



11 

^-^exp(elog(4fc„^6)) 



8k 



rid 



r ( - + 2zE^,ndrG 



log (4fc„^6) 



(A13) 



where we used the property r(e/2) = 2r(l + e/2)/e. Finally, by exploiting 
Eq. (]A7[) we have 



/oo 
dz\Ky2+2zEr,,„^rGi4:kn^z)\ 



TT 



rlog(4A;„» (A14) 

8fc„,cosh(2E„,„,rG) ^^ "'^ ^ ^ ^ 

and, by substituting Eqs. (IXSl) . (lAlOl) and (I A 141) in Eq. (IX9|) we obtain the 
following expression of Nn,' 



(» 



nd 



AtO) = A^O) 

n,nd Ud 



An 



rG[l-81og(4A;„,6)]' 



(A15) 



Since, at the end of the calculations the limit 6^0 will be performed, we 
can give the expression of Nn^ in this limit: 



rid 



n 



2rG log (kn^b) 



knj) ^ 0. 



(A16) 



In the same limit an easy expression of the density of the energy levels 
Q {En^rid) can be obtained. In fact, this quantity is defined as 



Q iEn,nd) = 



dn 



dE. 



n,nd 



Now, if kn.b — > then Eq. flA6p becomes simply 



IT 



En,nd'r'G log {knJ)) = Tl- 



knj) —> 



(A17) 



(A18) 



27 



and the density of the energy levels does not depend on the energy itself: 

2rG 

g {En,na) = Qn^ = log (^n^^) K^b ^ 0. (A19) 

TT 

Finally, with this definition the normalization factor A^„^ in the limit fc„^6 — > 
can be written as 

iV(>) = -^ K^b -. 0. (A20) 
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